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Abstract 

We introduce a rotation invariant short distance cut-off in the theory of an 
ideal fluid in three space dimensions, by requiring momenta to take values in 
a sphere. This leads to an algebra of functions in position space that is non- 
commutative. Nevertheless it is possible to find appropriate analogues of the 
Euler equations of an ideal fluid. The system still has a hamiltonian structure. 
It is hoped that this will be useful in the study of possible singularities in the 
evolution of Euler ( or Navier-Stokes )equations in three dimensions. 

O : 

1 Introduction 

| A square-integrable function in Euclidean space can be represented as a Fourier 

Oh 1 transform 
Jl f(x) = / f(k)e 2 ^dk 

a The wavenumber k also takes values in Euclidean space, which is naturally 
thought of as the dual of the original space. If the function is smooth its Fourier 
transform will decay faster than any power in the dual space. 
£N1 ■ If we consider instead functions on a lattice / : Z d — > C, the Fourier trans- 

form is a function on the torus / : T d — > C . In the language of solid state physics, 
^} • the fundamental domain of this torus is the 'Brillouin zone' of momenta of an 

electron in a periodic potential. Because there is a smallest possible length ( 

■ the distance between nearest neighbors in the lattice) there is a largest possible 
| momentum (half the diameter of the Brillouin zone). 

■ Thus there is a reciprocal relation between the smallest distance allowed in 
| space and the largest allowed wavenumber. It is analogous to the uncertainty 

principle of quantum mechanics. Indeed, it is the uncertainty principle, once 
it is accepted that particles are represented by waves. 

When studying a partial differential equation it is often useful to impose such 
a smallest possible length scale, at the cost of introducing some non-locality in 
the problem. This happens if we replace space by a lattice to discretize the 
PDE to solve it numerically. As noted above, this imposes a limit on the mag- 
nitude of the largest possible wavenumber. In quantum field theory, such a 
cut-off in momentum ('regularization') is needed to control the divergences that 
arise inevitably. After this regularization, we would study the dependence of 
the parameters ( coupling constants) on the cutoff implied by the condition 
that physical observables be independent of the cutoff. This is the still mys- 
terious method of 'renormalization'. Applications of quantum field theory to 
high energy physics as well as critical phenomena require cutoff procedures that 
preserve the symmetries of the original field theory. 

The lattice method is not always the best regularization, as it breaks rotation 
invariance. Also, replacing space by a lattice of points introduces a lack of 
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smoothness of functions. Consequently, numerical methods of solving PDEs 
suffer from spurious instabilities. 

Is it possible to introduce a smallest possible length in space without break- 
ing rotation invariance, and while maintaining smoothness of functions in space? 
Some proposals of this kind have appeared in theories of quantum gravity and 
in particle physics [TJ HI I3j (where the symmetry of interest is Lorentz invariance 
rather than rotation invariance). The price we pay for this is a kind of fuzzy- 
ness in space, where its co-ordinates become non-commutative. It is no longer 
possible to locate particles with infinite precision; the algebra of functions on 
space is replaced with a non-commutative algebra. These techniques could be 
useful in the theory of PDEs and QFT; and also of practical use in solving PDEs 
numerically. 

In this paper we will investigate the consequences of introducing such a small- 
est possible length scale into fluid mechanics in three dimensions. There is some 
earlier work on two dimensional fluid mechanics [U El O Although there 
are some situations where it is possible to limit fluid flow to two dimensions, 
the vast majority of phenomena of interest are in three dimensions. Some of 
the deepest problems in physics ( turbulence) and in mathematics (existence of 
solutions to Navier-Stokes equations) are in three dimensional fluid mechanics. 
A fundamental phenomenon is that information flows into large structures in 
the fluid from small distance scales (causing apparent randomness of the large 
scale degrees of freedom) while energy flows into small scales (dissipation due to 
turbulence and viscosity). Mathematically three dimensional Navier-Stokes 
is a 'hard' PDE because it is 'supercritical': the nonlinearities become stronger 
at small distance scales, making it impossible to know using present techniques 
whether solutions remain smooth for all time. Thus it is crucial to understand 
the scale dependence of fluid mechanics. 

Experience from quantum field theory suggests that we must first replace the 
Navier-Stokes equations with a 'regularized' version, in which there is a short 
distance cutoff. It should be possible to understand the regularity of this cutoff 
version. As pointed out by Tao, the problem of establishing is now shifted to 
studying the limit as this cutoff goes to zero. The experience of quantum field 
theory suggests that a cutoff that preserves rotation invariance is needed to have 
a 'renormalization' theory of the limit as the cutoff is taken to zero. At least it 
is a worthwhile avenue of study. 

While removing the cutoff is a great mathematical challenge, the cutoff the- 
ory itself could of some interest in physics. After all, the equations of fluid me- 
chanics are an approximation valid for an average of a large number molecules. 
A 'fuzzy' version of fluid mechanics would describe even larger scale motion, 
which averages over fluid elements themselves. Such a 'mesoscopic' theory may 
be what we need to understand many physical phenomena, such as the stability 
of large vortices. Because we averaged out fluid motion, the effective theory of 
large scale motion can be non-local; i.e., the we get integro-differential rather 
then differential equations. 

Computational Fluid Dynamics is important to many engineering applica- 
tions from weather prediction to the design of aircraft. Typically [TD], space is 
divided into a finite number of cells. The PDEs are turned into finite differ- 
ence equations that arc solved numerically If the size of the cell can be made 
small enough this can give a good approximation to the real flow. However, 
the number of cells is limited by the memory of the computer. If the region 
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of space is large the size of a cell can be too large. For example in weather 
prediction, a cell is several kilometers in size. This means not only that you 
miss phenomena within such cells, but also that predictions are limited in time. 
Given enough time the small scale will affect the large scale motion. In the case 
of the atmosphere the limit is about ten days beyond which predictions of the 
weather become unreliable with even the largest computers. 

Thus, a method that imposes a smallest possible length, and a largest possi- 
ble wavenumber, without breaking symmetries could help us in mathematical, 
physical and engineering approaches to fluid mechanics. 

Numerical simulations of quantum held theory [3] by the non-commutative 
regularization sometimes runs into trouble, because of a kind of anomaly. A 
remnant of non-commutativity survives even in the limit as the regularizaton 
is removed. Thus, the limiting theory is not what one would have hoped to 
have originally, ft can be asked [llj whether such problems can arise in classical 
held theories such as fluid mechanics. At least in two dimensions, the answer 
is no. It has been rigorously established [8] that the limit of the solution of the 
regularized theory ( on a torus) is indeed the solution of Euler's equation. For a 
more detailed and updated discussion of these issues we refer the reader to the 
author's webpage [TT] . 



2 The Euler Equation in Three Dimensions 

We now summarize the relevant facts about fluid flow in three dimensional fluid 
flow [H[T3]. 

The Euler equations of an incompressible inviscid fluid are 
d 

—v + (v ■ V)v = -Vp, div v = 0. (1) 

We can eliminate pressure p by taking a curl. Defining the vorticity uj = curl v 
and using (v ■ V)i! = u> x v + ^Vv 2 , we get, 

d 

—uj + curl [cj X v] = 0. (2) 

We can regard u) as the basic dynamic variable [Hj of the system, since v deter- 
mined by it uniquely as the solution of the equations 

div v = 0, curl v = uj. (3) 



2.1 Clebsch Variables 

It was noticed by Clebsch that any solution to div uj = can be written as 

uj = VA x (4) 

for a pair of functions A, [i : R? — > R. The Euler equations take a simple form 
in these variables [H]. They are constant along stream lines: 

^ + v ■ VA = 0, ^ + v ■ Vm = 0. (5) 

1 We assume that the velocity vanishes at infinity; this is necessary to have finite energy. 
Or, we could study fluid flow in some domain in R 3 with appropriate boundary conditions. 
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Here we are to regard velocity as a function of the Clebsch variables: 

u = P(AVju). (6) 
Also, P is the projection of a vector to its divergence free part: 

P(u) =u + V4>, V • u + V 2 = 0. (7) 
This follows from the identity 

VAxV/j = Vx [AV/i] = V x v. (8) 

2.2 Some Technical Remarks 

1. We will lose some global information about the flow in this parametriza- 
tion; for example the 'kinetic helicity' (a conserved quantity) will vanish if 
A, \i are globally defined: j w X vdx — 0. To avoid this, we must let (A, //) 
be co-ordinates on a two dimensional manifold more general than R 2 . But 
we postpone such issues for now. 

2. There is a 'gauge invariance' to this parametrization: under canonical 
transformations in the pair of co-ordinates (A,/x), vorticity is unchanged. 

3. J. Marsden and A. Wienstein.15] have given a nice geometrical interpre- 
tation of these variables. The space of pairs of functions (A, /z) is a sym- 
plectic vector space using the contraction < A,/i >= J A/id 3 x given by 
the volume form. The group of incompressible diffeomorphisms act on 
this symplectically. The Clebsch bracket co = [[A, fi]] := VA x V/i is the 
moment map of this action. This is anlogous to the formula L = r x p for 
angular momentum. 

4. We propose another interpretation that suggests a natural deformation. 
Let A, [i : G* — > R be pair of functions on the dual of a Lie algebra G. Then 
dX and d/i can be thought of as functions : G* — > G. Thus it makes sense 
to define the commutator [dX, dfj]. This is a map V <S> V — > V <E> G, where 

V is the space of real functions on G. As a vector space we can identify 

V = U(G), the universal envelope of G. In the present case G = SU(2) 
and its Lie algebra dual is R 3 . In a later publication we will discuss the 
deformation of this universal envelope into a non-co-commutative Hopf 
algebra ( " quantum group" ) which provides a natural regularization of 
both short-distance and large distance effects. 

2.3 The Hamiltonian in the Clebsch Parametrization 

Another advantage of the Clebsch variables is that ideal fluid flow can be thought 
of a hamiltonian system in which they are the canonical co-ordinates: 

{X(x),fi(y)} = S(x,y), {X(x), X(y)} = = {^(x),fx(y)}. (9) 

The hamiltonian is just the total kinetic energy of the fluid: 

H =^<v,v>=^< AV/i|P|AV/i > . (10) 
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It is not difficult to verify that the Eulcr equations in Clebsch form follow from 
this hamiltonian and Poisson bracket. 

It will be useful to write the hamiltonian in terms of Fourier transforms, 



\(x) = / \(k)e 



2-nik-x 



dk, 



X*(k) = A(-fc). 



(11) 



We get 



{A*(A0,£(p)} 



S(k,p) 



(12) 



all other Poisson brackets being zero. Also, 



H =2j dkd P dk ' d P s ( k +p-k'-p') 
A*(fc')M*(p')A(fc)Ab)(27r) 2 \ P '-p- 



p' ■ (k + p)(k + p) ■ p 
\k+p\ 2 



3 Bounded Momentum Space 

A way of getting a configuration space of finite extend for a particle is to imag- 
ine that it is moving on a 3-sphere. This preserves rotation invariance, unlike 
periodic boundary conditions that put the system in a torus. In the limit of 
infinite radius we will get back Euclidean space. Its phase space will be T*S 3 ; 
using the group identification S 3 w SU(2) , we get T*S 3 w S 3 x R 3 . 

How will be put a cut-off on largest possible wavenumber? Does it make 
sense that momentum space is 5 3 ? 

Any function / : R 3 — > C has the Fourier representation 



The product of two functions corresponds to the convolution defined through 
the addition of momenta: 



If we change the composition law for momenta the rule for multiplication of 
functions will change as well. If momenta are valued in the non-abelian group 
SU(2), the multiplication of functions in space becomes non-commutative. 

Suppose we are given a smooth, locally invertible, function k : SU(2) — ► R 3 
such that k(l) = 0. We can then get a cutoff version of the Fourier transform 




(13) 




(14) 




(15) 



Here dg is the invariant (Haar) measure on SU(2) normalized so that 




(16) 



The inverse of this is 




(17) 
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where J(g) = det || is the Jacobian of the transformation k — > g. 

We also define a *-product on functions on R 3 , induced by the group multi- 
plication on 5/7(2): 

e 27rifc(gi)-:r ^ e 2ixik{g 2 ) -x _ e 2^ik(g 1 g 2 )-x 

More explicitly, 

fi * / 2 (x) = [ h{y)h{z)K{x,y,z)dydz (19) 



where 



Jf(a;,y,«) = ^ d gi dg 2 e 2lTt ^ 9l92 > x ' k ^- v ~ k< - S2 >^ 



(20) 



Associativity of this multiplication follows from that of group multiplication on 
SU(2). But it is not commutative. 

The momentum operators are, in this picture, just multiplication: 

PJ(k) = kj(k). (21) 

Hence 

[^,^•1=0. (22) 

Ordinarily the position operators would be the differentiation with respect to 
k. Instead, they will be the generators of the group action: 

A (23) 

so that 

[Pi,Xi]=i${P) (24) 

The components of this vector field on the group are determined in terms of the 
Maurer-Cartan forms: 

e j {k)uj{k) = Si, g- x dg = w\dk^. (25) 

It follows that £,j{k) = 5 % k + O(fc). Thus for momenta small compared to the 
cutoff, the position and momentum operators satisfy the usual canonical com- 
mutation relations. 

We can see that The position operators no longer commute; instead they 
satisfy the commutation relations of the Lie algebra. 

[X\X j ] =ie%X k . (26) 

Therefore it will not be possible to locate the position of a particle with infi- 
nite precision. This fuzziness in position means that in effect we are averaging 
physical quantities over a small region whose size a is a kind of cutoff in space, 
of the order of the inverse of the diameter of the momentum space. 
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An example of such a function k is (the inverse of ) the stereographic co- 
ordinate system 

1 — a 2 k 2 2iaa ■ k 
9= l + a 2 k 2 + l + a 2 k 2 ' ( ' 

where a is a constant with the dimensions of length. In these co-ordinates, the 
point at infinity in the variable k corresponds to g = — 1. Thus we must require 
all functions to tend to a constant as \k\ — > oo. The formula for the modified 
Fourier transform can be written as 



f{x) = J /»< 



2-nik-x 



a 3 dk 
(1 + a 2 fc 2 ) 3 



Pig) 



arccos 



3 jr 

since dg = ^Ji^-iyi in these co-ordinates. 
If we set 

we can solve for k to get 

\k\ 



tig 



1 1 / >. , i tan \p{9) . 

-tan-p(fif), k = -- . \ tvga 

a 2 2a sin p[g) 



(28) 



(29) 



(30) 



In the limit as a — > the cutoff Fourier transform and the ^-product reduces to 
the usual ones. To be specific we will assume from now on that we are using 
the stereographic parametrization. Up to second order terms in a, 



g ~ 1 + 2aia ■ k, k{g) w — — tvga, fc(.9i<?2) ~ k\ + ki + aki x + ■ ■ ■ (31) 

Other choices of the function k will give fluid equations that differ by terms 
higher order in a. Such an ambiguity in the choice of regularization is common. 



4 Euler's Equation With Short Distance Cutoff 



It is no problem generalizing the Poisson brackets; they just involve delta func- 
tions. 



{A*(5i),A(ff 2 )} = %r 1 ff 2 ) 



(32) 



The other brackets {X(gi), A (52)} and fi(gi), ^(52)} remain zero. 

The hamiltonian needs some work. We should replace k,p,k',p' by elements 
of SU{2). It is useful to denote 



1 

4^2 



< 5i,52 >= 

As a -> 0, < gi, g 2 fci • k 2 . 
This way we can get a formula 

H = \ J dgidg2dg '^ dg ^ 5 ^ g '^ 1 9 ^ 



tr 3i.92 - 2 



(33) 



7 



A*(<?i)A*(<7 2 )A(<7i)/i( ff2 ) <g' 2 ,g 2 > 



< 92,9192 >< gW2,92 > 

< .9132,5152 > 



Thus, in principle we have a cutoff version of three dimensional fluid me- 
chanics. Hamilton's equations follow: 



It would be interesting to see if these equations have a smooth time evolution 
given smooth initial data. 

This is not quite ready for a numerical applications, because we still have an 
infinite number of variables. We have a cutoff in momentum space but position 
space is still infinitely large. A cutoff in both spaces is necessary to get a finite 
number of degrees of freedom. This can be done, but the resulting algebra of 
functions is not any more related to a group. It is still associative and is related 
to a quantum group ( Hopf Algebra). We hope to return to these points. 
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